Josephson lattice structure in mesoscopic intrinsic Josephson junctions by means of 

flux-flow resistance in Bi 2 Sr 2 CaCu208 +( 5 
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Dynamical nature of the Josephson vortex (JV) system in Bi2Sr2CaCu20g+5 (Bi2212) has been 
investigated in the presence of the c-axis current with magnetic field alignments very close to the 
ab-plane. As a function of magnetic fields, the c-axis JV flux flow resistance oscillates periodically 
in accordance with the proposed JV triangular structure. We observe that this oscillating period 
becomes doubled above a certain field, indicating the structure transition from triangle to square 
structure. This transition field becomes lower in junctions with smaller width perpendicular to the 
external field. We interpret that this phenomena as the effect of the edge deformation of the JV 
lattice due to surface current of intrinsic Josephson junctions as pointed by Koshelev. 
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I. INTRODUCTION 

Highly anisotropic high-T c superconductors such as 
Bi 2 Sr 2 CaCu 2 8+< 5 are well described as stacks of weakly 
coupled Josephson junctions^, where charge coupling be- 
tween Cu0 2 double-layers is not negligible^. This charac- 
teristic feature is responsible for properties of so-called in- 
trinsic Josephson junctions (IJJ) like multiple branches in 
current-voltage characteristics and longitudinal Joseph- 
son plasma modei Appearance of various vortex phases 
under the c axis field indeed manifests the weak coupling 
between the layers. 

In the magnetic field parallel to the ab plane (_L c axis), 
the situation becomes different because of inductive cou- 
pling between Josephson vortices generated by the paral- 
lel magnetic field^. In high fields above > 1 T, Josephson 
vortices occupy all block layers and form the Josephson 
vortex lattice which demonstrates the strong interaction 
between vortex arrays in the neighboring layer* 5 . This 
is the one of the most striking difference of Josephson 
vortices in IJJ from ones in isolated single Josephson 
junctions. Furthermore, many characteristic features of 
Josephson vortex state in parallel magnetic fields have 
been found in the I — V characteristics^, the Josephson 
plasma resonance^, and the vortex phase diagram^. 

Recently, Ooi et a! 9 found periodic oscillation of the 
Josephson vortex flow resistance in small Bi2212 single 
crystal as a function of magnetic field. They claimed oc- 
currence of the coherent flow of the triangular Josephson 
vortex lattice. Subsequently, Koshele-\ii£ and Machida 11 
have reproduced the oscillation by calculating the flux 
flow resistance by taking the surface barrier effect into 
consideration. This means that the surface barrier is not 
so strong to destroy the bulk properties of the Josephson 
vortices, where it prefers to form the triangular lattice 
even in the vicinity of sample edges. This is in contrast 
to the vortex state in fields parallel to the c axis where 
vortex arrangements in the vicinity of sample edges are 
dominated by the surface barrier effect. 

In this paper, we present experimental results of the 
Josephson vortex flow resistance in Bi 2 Sr 2 CaCu 2 08+,s 
micro-structured crystals with dimensions of 1—10 /xm. 



In the smaller samples and at higher fields, we have found 
periodic oscillation of the Josephson vortex flow resis- 
tance. This indicates that the square vortex lattice is 
realized. This experimental observation is attributed to 
the effect of the surface barrier, which becomes dominant 
in smaller samples and overcomes the bulk properties in 
the short stack regime. 



II. RESULTS AND DISCUSSIONS 

We made in-line symmetric junctions^ to avoid vari- 
ous contact problems accompanied by the transport mea- 
surements of intrinsic junctions. Bi 2 Sr 2 CaCu 2 0g+,5 sin- 
gle crystals were grown by the travelling solvent floating 
zone method. A cleaved thin (~ 10 fim) crystal was care- 
fully cut into narrow strips with a width of ~ 100 /jm by 
a sharp knife. The narrow strip was put on a MgO sub- 
strate with four gold electrodes which were patterned by 
a vacuum deposition with a mask and the center of the 
strip was shaped into mesoscopic junctions as shown in 
the inset of Fig. njwith the focused ion beam (FIB) ma- 
chine SMI2050MS (SII NanoTechnology Inc.). We made 
junctions from eight different crystals from A to H and 
FIB processes were repeated three times in crystal H af- 
ter measurements, resulting in eleven samples in total 
with different dimensions as listed in Table [I] 

By applying the c axis current between two current 
electrodes, the current is concentrated into small junc- 
tion area only. The voltage between the voltage elec- 
trodes comes from only the small junction area because 
the other parts of the sample still remain in the supercon- 
ducting state. The temperature dependence of resistivity 
of a junction is shown in Fig. ^ in which a clear super- 
conducting transition and zero resistance below T c are 
obtained. External magnetic fields were applied by the 
split pair superconducting magnet which can generates 
horizontal magnetic field up to 80 kOe. Rotating mag- 
netic field around the ab-plane with a precise rotator, the 
c-axis resistance suddenly increases in the vicinity of the 
a6-plane as shown in the inset of Fig. ^ This resistance 
is attributed to the Josephson vortex flow resistance and 
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TABLE I: List of samples which we measured. Numbers fol- 
lowing capital letters refer to sizes of junctions perpendicular 
to the magnetic field and the c-axis. 





w [/im] 


1 [/xm] 


t [fim] 


T c \K] 

c L J 


A94 


9.4 


10.3 


0.3 


86.6 


B19 


1.9 


11.7 


0.5 


90.0 


C22 


2.2 


8.9 


0.14 


83.6 


D49 


4.9 


17.8 


0.08 


88.2 


E44 


4.4 


9.5 


0.15 


84.7 


F18 


1.8 


9.9 


1.14 


87.2 


G20 


2.0 


10.2 


0.9 


84.3 


H41 


4.1 


5.5 


2.0 


86.9 


H55 


5.5 


4.1 


2.0 


86.9 


H38 


3.8 


5.5 


2.0 


86.9 


H30 


3.0 


5.5 


2.0 


86.9 


H23 


2.3 


5.5 


2.0 


86.9 



the lock-in state without pancake vortices is realized in 
this angle range. We can set the field alignment parallel 
to the afr-plane with an accuracy of 0.01 degree. 

Figure |2] (a) represents field dependence of the c-axis 
resistance R c at 60 K in H55. External field was applied 
parallel to the a6-plane and swept up to 50 kOe with a 
constant c axis current (1 /iA) at a constant temperature. 
With increasing field from zero, the resistance smoothly 
increases up to 8 kOe, then begins to oscillate with oscil- 
lation center being monotonically increasing. The oscilla- 
tion period H p is constant up to 30 kOe, H p = 1.25 kOe, 
and the amplitude is maximized around 25 kOe, above 
which the oscillation gradually diminishes with increas- 
ing field. However above 35 kOe, the other oscillation in 
resistivity with different H p develops and the oscillation 
of Hp = 1.25 kOe seems to disappear above 40 kOe. 

H p below 30 kOe is very close to the field, which cor- 
responds to add a half flux quantum 4>o/2 to a JV array 
in a Bi2212 block layer as Ho/2 = (j) /2ws = 1.25 kOe 
for H55, where s — 15A is the periodicity of CuC>2 dou- 
ble layers of Bi2212. This oscillation is quite similar to 
the previous results observed by Ooi et al£, where they 
argued that this is due to surface barrier effect for JVs 
forming the triangular lattice. The surface barrier for the 
triangular JV lattice is considered to be enhanced when 
the number of JVs in a layer (H/Hq) corresponds either 
an integer or a half odd integer because the periodicity of 
the JV lattice matches with the junction width for both 
two cases. These two matching result in the oscillation 
having period of Hq/2 and the bottom at H/Hq being 
half integers. 

Hp above 40 kOe was found to be 2.49 kOe, which is 
twice as the period below 30 kOe and very close to the 
Hq value of 2.50 kOe. With increasing field in transient 
region from the Hq/2 oscillation to the Ho oscillation, the 
amplitude of Hq/2 oscillation diminishes and the ampli- 
tude of the Hq oscillation develops, resulting in the oscil- 
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FIG. 1: Temperature dependence of the c-axis resistivity at 
H = in E44. It is found that the crystal is slightly over- 
doped . The lower inset shows angular dependence of the 
c-axis resistivity of J = 50 A/cm 2 at H = 1 T in the vicinity 
of the ab plane. Arrows indicate directions of change in angle 
and slight hysteresis was found. Upper inset depicts SEM 
image of the sample. A thick bar corresponds to 10 /im in the 
image. 

lation amplitude is larger than one of the Ho/2 oscilla- 
tion. The Ho oscillation starts at lower fields in samples 
with smaller width as shown in Fig. (b) where Ho /2 
oscillation is hardly seen even at the lowest field, while in 
the largest sample A94 H /2 oscillation was observed up 
to the maximum field range. We now understand that 
the Hq/2 oscillating behavior observed by the previous 
measurement by OOi et al. is due to the large sample 
width of more than 20 microns^. 

Critical current density J c as a function of H» of H41 
is shown in Fig. |3 which is extracted from I — V char- 
acteristics at various H\\ with a criterion of V = 1 mV. 
One finds that J C (H) below 35 kOe oscillates with pe- 
riod of Hq/2 and has a local maxima at H/Ho being half 
integers, whereas J C (H) above 42 kOe shows oscillation 
with the period being Hq and local maxima at H/ Hq be- 
ing half odd integers. Magnetic field dependence of J c 
similar to our results in low fields has been expected the- 
oretically by taking boundary deformation of JV lattice 
into consideration^. Also, J C (H) in high fields is similar 
to magnetic field dependence of J c in single Josephson 
junctions in a sense that the points where the magnetic 
flux inside the junction equals to the half odd integer 
flux quanta corresponds to the local maxima of the crit- 
ical current. The Fraunhofer-like oscillation of J C (H) is 
interpreted that the Josephson vortices form a square 
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FIG. 2: Field dependence of the c axis resistivity in H55 (a) 
and B19 (b) and threshold field Hth of various samples (c). 
Arrows in (a) and (b) point respective Hth and length of error 
bars in (c) corresponds to the period of the Ho/2 oscillation 
of the sample. 



lattice where all layers are identical to edges. 

Such a Fraunhofer oscillation in J C (H) of multilayered 
system like Bi2212 has been expected when the condition 
w -C A j = 7s is satisfied 13 . However, Aj is estimated to 
be in a range from 0.075 to 1.5 /im for slightly over-doped 
Bi2212, which is much smaller than w of this sample. 
This discrepancy seems to occur due to inhomogeneity 
of the current, which is more realistic in the real speci- 
mens of high-T c superconductors. It is known that the 
Bean-Livingstone surface barrier is pronounced in high- 
Tc superconductors, so that vortices prefer to stay at the 
edge of the sample. This boundary effect can also explain 
Ho/2 oscillation in J C (H) as discussed by Koshele\ii£. 

The boundary effect for the JV lattice should be more 
pronounced in smaller junctions because the number of 
vortices at both ends of vortex arrays inside a layer (usu- 
ally two) , which are most seriously affected by the bound- 
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3: Field dependence of J c extracted from a measurement 
- V characteristics. 



ary effect in smaller junctions has a greater ratio to the 
number of total vortices of the array than in the case of 
larger samples at constant fields. Consider a Josephson 
junction with two JVs as an extreme case. In p — H curve 
as shown in Figs. |U(a) and (b), we define threshold field 
H t h above which local minima in the vicinity of H / ' H a 
being an integer cannot be resolved (indicated by arrows 
in the figure). It was found that H t h is lower in samples 
with smaller w from p—H curves in samples with various 
w as plotted in Fig. 0(c). These results suggest that the 
square JV lattice is more favorable for smaller junctions 
and higher fields. 

Now, let us consider how vortices inside the junctions 
interact with edges and each other in transition from tri- 
angular lattice to square lattice. At low fields (although 
high enough to form dense lattice), the JV system is 
dominated by the interaction between vortices in adja- 
cent layers because the separation between JVs inside a 
layer is too far to convey the boundary deformation to 
all vortices in the layer. With increasing magnetic field, 
vortices inside the layer become closer and closer while 
the separation between layers does not change and the 
interaction between adjacent vortices in a vortex array 
and the boundary effect becomes stronger and stronger 
because JVs are more closely packed in the layer. The 
JV system is finally dominated by the boundary effect 
and turns to the square lattice in high field region at the 
cost of inductive coupling between layers. The inductive 
coupling is weaker in a smaller junction because the to- 
tal number of vortices in a layer is fewer, resulting in the 
transition from triangle to square at lower fields. 
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According to theoretical calculation by Koshelev, JVs 
with their center being closer to both edges than the 
lattice deformation length Ib — r ys(njs 2 B/\ // 2(f>o) devi- 
ate from positions extrapolated from the bulk region as- 
suming an undcformcd lattice. Here 7 and B are the 
anisotropy parameter and induced flux density in the 
sample, respectively. This means Josephson vortices in- 
side samples with w < 21b may not form the triangular 
lattice, so that the threshold field above which the sur- 
face barrier effect dominates the system is written as a 
function of w, 



Hi 



w (f>p 

7s 27T7S 2 



(1) 



A thick line in Fig. [21 (c) corresponds to Eq. Q with 
7 = 110. The agreement is rather good and the value 
of the obtained 7 is reasonable for slightly over-doped 
Bi2212. 



III. CONCLUSION 

We have investigated Josephson vortex lattice struc- 
ture in mesoscopic (~ /im) intrinsic junctions with the 



probe of the JV flow resistance along the c axis. For the 
small (w < 6 /zm) mesas, the JV lattice is found to change 
the structure from triangular to square lattice with in- 
creasing fields in low current limits. Transition field from 
the triangular to square lattice monotonically increases 
with w. This is attributed to the interplay between the 
surface barrier effect and the inductive coupling between 
Josephson vortex arrays in intrinsic Josephson junctions. 
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